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ABSTRACT 

In  this  paper  a delayed  SIR  model  with  exponential  demographic  structure  and  the  nonlinear  incidence  rate  is 
formulated.  We  show  if  the  basic  reproductive  number,  denoted,  R0 , is  less  than  unity,  disease  free  equilibrium  is  stable. 

Moreover,  we  prove  that  R0  > 1,  the  endemic  equilibrium  is  locally  stable  without  delay  and  the  endemic  equilibrium  is 

stable  if  the  delay  is  under  some  condition.  Finally  a numerical  example  is  also  included  to  illustrate  the  effectiveness  of 
the  proposed  model. 
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Bifurcation 

AMS  Subject  Classifications:  34C07,  34D23,  93A30,  93D20 

1.  INTRODUCTION 


In  recent  years,  more  and  more  delayed  models  have  been  investigated  during  the  study  of  epidemic  models 
[3,  5,  6,  8,  10,  11,  13,  14].  Hethcote  and  Van  den  Driessche  [4]  have  considered  an  SIS  epidemic  model  with  constant  time 
delay  which  accounts  for  duration  of  infectiousness.  Beretta  et  al.  [1]  have  studied  global  stability  in  an  SIR  epidemic 
model  with  distributed  delay  that  describes  the  time  which  it  takes  for  an  individual  to  lose  infectiousness.  Kaddar  et  al.  [7] 


considered  a delayed  SIR  model  with  a saturated  incidence  rate 


pmm 

1 + (X^S  ( t ) + OCrj.  (t) 


as  follows: 


*=A-W) 

dt  1 + CfjSV)  + 

dl__  0S(t-T)I(t-T) 

dt  1 + alS(t-T)  + a2I(t-T) 
dR 


~(ju  + a+y)I  (t), 


dt 


yI(t)-juR(t). 


j3S(t)I(t) 

The  characteristic  of  this  model  is:  the  saturated  incidence  rate which  includes  the  three 

1 + (X^S ( t ) + CCrJ.  it) 

forms:  j8S(t)I(t)  ( if  OCx  —OC2  =0),  — — ^ ( if  OC2=  0)  and  — — — ( if  (X{=0)  saturated  with  the  susceptible 


1 + axS(t) 


l + OCryl  (/) 
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S(t)  and  infective  I(t)  individuals.  The  inclusion  of  time  delay  into  susceptible  S (t)  and  infective  I(t)  individuals  in 
incidence  rate,  only  on  the  second  equation,  because  susceptible  individuals  infected  at  time  t — T is  able  to  spread  the 
disease  at  time  t.  Rihan  et  al.  [9]  considered  a qualitative  analysis  of  delayed  SIR  epidemic  model  with  saturated  incidence 

rate as 


1 + <jS(t) 

dS 


dt 


= rS(t ) 


V 


5(0")  0S(t)I(t-T) 
K ) 1 + <t5(0 


dl  _pS(t)I(t-T) 
dt  1 -\-(rS(f) 
dR 


-al(i)-al(t\ 


dt 


= aI(t)-SR(t). 


Here  parameters  r is  the  logistic  growth  rate,  K is  carrying  capacity,  <J  is  the  saturation  factor  that  measures  the 

1 * 
inhibitory  effect,  — is  the  incubation  period  [2,12],  CC  is  the  recovery  rate,  O is  the  natural  death  rate  due  to  causes 

T 

unrealed  to  infection  and  a is  the  infected  hosts  die  rate  which  includes  both  the  natural  death  rate  plus  the  disease  induced 
death  rate. 


In  this  paper  we  consider  a delayed  SIR  model  with  the  nonliner  incidence  rate 

analyze  the  stability  and  the  existence  of  Hopf  bifurcation. 

2.  DELAYED  SIR  EPIDEMIC  MODEL 


/3Sit-r)l{t-T)e 
\ + aS(t-z) 


-jUT 


We  also 


In  this  section,  we  consider  the  following  SIR  model  with  the  non  linear  incidence  rate 


j3S(t-T)I(t-T)e 
l + aS(t-z) 


Let5(?)is  the  number  of  susceptible  individual,  I(t ) is  the  number  of  infective  individual,  and  R(t)  is  the  number  of 
recovered  individuals,  then  we  have  the  following  model 


dS  cro 

— = rS(t) 

dt 


1- 


5(f)")  pS{t-T)I(t-T)e 


fIT 


dl  _ /3S(t  - t)1  (t  - r)e“'ur 


dt 


K J l + aS(t-r) 

{fi  + p+S)I(t), 


--//5(f), 


l + «5(f-  r) 


-j-  = pl(t)  - ^R(t). 

dt 


(2.1) 


The  parameter  r is  the  logistic  growth  rate,  K is  the  carrying  capacity,  jl  > 0 is  the  rate  of  natural  death  such 


that  r > jd  , S > 0 is  the  rate  of  disease  related  death,  p > 0 is  the  rate  of  recovery, 
the  parameters  that  measures  infections  with  the  inhibitory  effect. 


1 

— is  the  incubation  period  and  CC  is 
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The  first  two  equation  in  system  (2.1)  do  not  depend  on  the  third  equation,  and  therefore  this  equation  can  be 
omitted  without  lose  of  generality.  Hence  system  (2.1)  can  be  rewritten  as 


dS 

— = rS(t) 
dt 


1— 


Sit)  \ pSit-r)I{t-T)e 


K ) l + a.S(f-r) 
dl  _ pS(t-r)I(t-r)e-MI 
dt  l+tz^t-r) 


-fiS(t), 


-ift  + p+S)I(t). 


(2.2) 


Propostion:  For  the  model  system  (2.2),  there  always  exit  infection  free  equilibrium 


£0=(0,0),  El  = 


( J-p)K 


,0 


.If 


K(r  - fi)\_J3e  '“r  - ct(fi  + p+ 5)] 
r(/t  + p+S) 


(2.3) 


There  also  exits  an  endemic  equilibrium  E — (5  , / ),  where 


S =- 


ifi  + p + S) 


,1  = 


rS 


[/VU- 


f5c  — ocip  + p + S)  K ip  + p + S) 

3.  LINEAR  STABILITY  ANALYSIS 

The  characteristic  equation  for  the  model  (2.2),  is  given  by 


. . 2 rS  Ble t“+'°r  , 

ir-u) — r-A 

K (1  + aS)2 

pie -^X)t 

i 1 + aS)2 


pSe 


-(P+a)t 


1 + ccS 


j3Se 


-G u+a)t 


1 + aS 


(Z7  + p+S)  - A 


(3.1) 


Theorem  3.1:  E()  is  always  a saddle  point  and  there  can  not  be  total  extinction  of  the  system  (2.2). 
Proof  Using  (3.1),  the  characteristic  equation  at  E()  = (0,0)  reduces  to 

[Z-(r- p)][A  + (p+p+S)\  = 0 

Obviously  (3.2)  has  a positive  root  — r — Thus  EQ  is  always  unstable  (saddle  point). 


(3.2) 


Theorem  3.2:  The  infection  free  equilibrium  Ex  — 
unstable  when  R0  > 1,  and  linearly  neutrally  stable  if  RQ  = 1 . 


K(r-jU) 


,0 


is  asymptotically  stable  when  ^ an<^ 
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Proof.  Using  (3.1),  the  characteristic  equation  at  Ex  — 

j3K(r-jU)e-<Jl+A)T 


K(r-jU) 


V r 


,0 


IS 


(r-jU)-A 


0 


1 + oc 


K(r-ju) 


l 3K(r-ju)e 


-(ju+A)t 


1 +a 


K(r-ju) 


— ( fi  + p + S)  — A 


= 0 


This  implies 


x+iM+p+S)li-rAl^^L 

^ H n r + aK{r-p) 


:0. 


(3.3) 


The  two  roots  of  (3.3)  are  real  and  negative  if  RQ  < l(whenT  = 0).  The  equilibrium  Ex  is  asymptotical  stable. 
WhenT  > 0,  we  suppose  (3.3)  has  a purely  imaginary  root  X — COi , then  separating  real  and  imaginary  parts,  we  have 


-<d2  +(r- //)(//  + p + S)  = {r-fi){fi  + p + S) 


rR^  + a(r  - p)K 
r + a(r-  p)K 


cos  COT 


(3.4) 


, +a(r-p)K  . 

+ co([i  + p + S) — — — sin  cor 

r + a(r-p)K 


[{fi  + p+  S)  + (r - fi)\co  = a{fi  + p+  coscot 

r + a{r-fi)K 

, \r  ri?o  +a(jr-fi)K  . 

-(r- //)(// + p+<5)  — 5 — srnrnr 

r + a(r  - p)K 


Hence 


(3.5) 


or  =(fi+p  + S)2 


rR^+ap-pKV  t 
r + a(r  — p)K  J 


(3.6) 


When  Rq  < 1,  then  there  are  no  positive  real  roots  CO. 
This  is  complete  the  proof. 


1 


Theorem  3.3:  (i)  The  endemic  equilibrium  E is  asymptotically  stable  if  1 < R,.  <2  H , when  T — 0. 

aS 


(ii)  When  T > 0,  R()  >1  suppose 


1 


1 + aS 


i--L 


1 


1 +aS 


i-i- 


2(jU  + p + 8 ) 
r- fi 


fi  2 ) 

< 

V Ro) 

, then 
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there  exit  Tx>  0 such  that  T G [0,  Tx  ] the  endemic  equilibrium  E is  stable,  and  unstable  when  T > Tv 
Proof  The  characteristic  equation  for  the  system  (2.2)  at  E = (S  , / ) is  given  by 


, , 2 rS*  J3l*e-(M+X)T  . 

( r-ju ) — 5^; — A 

A1  (1  + aS  )2 

(1  + arS*)2 

This  implies 


/?5V 


* -(/<+A)r 


1 + orS 


-(//+/l)T 


1 + aS 


■ + p + — A 


= 0. 


(r-//)- 


2 (r-ju)  r-ju 


R 


0 

r-// 


1 + aS 


1 + aS 


R-o  J 


?~Xr-A 


-At 


~(ju  + p + S)e  Xt 
(ju p 8^)  [e  — 1 j — A 


= 0. 


r+  -(r-M) 

+(/i  + p + S) 


1 — 

^0  J 


+(fi+p  + 8)(\-e-u)  + 


r~M 


1 + aS 


1-  — 

V 


\-2-' 

A 


(i-^)+T 


r-fi 


+ aS 


1-— 

V 


Introducing 


Qi  =(r-n ) 


1- 


V iV0  y 


Rn 


r - ji 


1 + aS 


1- 


V ^0  J 


Rn 


Then  the  characteristic  equation  (3.7)  can  be  rewritten  in  the  form 

/L  + + P + [l^|A  + ^ ]s  — 0. 

Where 


',-C 

v Rq  j 


Pi  = ~Q\  + Q3  = ~(r  ~ P) 

P2  = -QiQ3  = -O  - M)(M  + P + #) 


+ (u  + p + 8), 


v Roj 


qx  - —Q3  + Q2  - -(//  + p + 8)  + 


( r~P ) 


1 + aS 


^oy 


q2  ~ Q1Q3  Q2Q3  ~ 


i-A+.  i 


i?0  1 + aS 


^oy 


C3  = (//  + /?  + 8). 


(r-ju)(p  + p + 8). 


(3.7) 


(3.8) 


(3.9) 
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When  T = 0,  the  characteristic  equation  (3.9)  becomes  X2  + axA  + a2  = 0.  Where 


<h  = Pi+<h  =(r-fi) 


1 


. 2 1 

-1  + — + - 


R0  1 + aS 


V J 


, a2  = p2  + q2—  Q2 Qx  > 0.  Therefore,  if  c/  > 0 i.e. 


R0  < 2 H , a2  > 0 (when  R0  > 1 ),  then  by  the  Hurwitz  criterion,  we  can  know  the  endemic  equilibrium  E is  stable. 

aS 


WhenT  > 0,  we  suppose  equation  (3.9)  has  a purely  imaginary  root  A = CQi,  the  separating  real  and  imaginary 


parts. 


CO  - p2  = q2  COS  COT  + qxG)sm  COT , 

= g2  sin  cot  - qxcocos  cot . 
Hence 

to4  +a3fi>3  +a4  = 0 

Where 

a3  = Pi2 -2p2-ql2 
= (r-jU) 

a2  ~ Pi  —q22<® 


(3.10) 


(3.11) 


(3.12) 


1 

V R0J 


(r-jd)1 


1 + aS 


i-L 

v ^oy 


1 


1+aS 


Vi' 

v ^oy 


r~P 


;{P  + P + 3) 


If  a2  > 0,  i.e. 


1 


1 


1 + aS 


'i--C 

^oy 


2 (//  + p + <?) 

r — ju 


< 


/ \ 2 

i-2 


*oy 


, then  equation  (3.12) 


1 + tfS 

has  at  least  one  positive  root,  say  C0{  > 0. 

Now,  we  turn  to  the  bifurcation  analysis.  We  use  the  delay  T as  bifurcation  parameter,  let  A(t)  = Y(t)  + iC0(r) 
be  the  eigenvalue  of  equation  (3.9)  such  that  for  some  initial  value  of  the  bifurcation  parameter  Tx , we  have  y(Tx ) = 0 and 
CO(Tx ) = COv  From  (3.10)  and  (3.11)  we  have 


tx  = — arccos 


(ga-AgiW  ~P2<h 
9iW  +q{ 


2 in  . , _ 

h , j — 0, 1, 2. 


(3.13) 


Differentiate  w.r.t.  Tx , we  get 


dV 

dr  j 


2^  + Pi 


4. 


-A(A  +^,A  + />2)  tf2A  92(^A  + ^f2)  A 


(3.14) 
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Thus 


p_ 

q: 

Where 

P = q2COl4  +a\2(2q2)  + q2(p2  -2  p2)  + p2q2,Q-^p20)2+(0)2  - R2)2][<?22  + q2C02~^  which  are 
positive  when  RQ  > 1 . 


sign  < Re 


dX\ 


= sign  < Re 


A=coi 


dX 

dr 


=sign 


Thus  we  have  i Re i >0,  by  continuity  the  real  part  of  M T)  becomes  positive  when  T > T,  and  the 

dT  L_, 


steady  state  becomes  unstable.  A Hopf  bifurcation  occurs  when  T passes  through  the  critical  value  T, . 


4.  EXAMPLE 


In  this  section,  we  present  some  numerical  results  of  system  (2.2)  at  different  T of  supporting  the  theoretical 
analysis  in  section  3 

We  take  the  parameters  for  endemic  equilibrium  without  delay  as  follows: 

5°  =25,/°  = 4,R°  = 6, r = ,4,/3  = 0.2,a  = 0.09, p = 0.3, p = 0.3, S = 0.2,t  = 0,K  = 60.  We  have 

R0  = 2.4  > 1, 5 —6.25,1  =0.4557,  then  by  theorem  3.3  (i),  Sit)  and  7(f)  approach  to  their  steady-state  values 

without  delay,  the  disease  will  be  exist. 

We  give  the  parameters  for  endemic  equilibrium  with  time  delay  in  system  (2.2)  as  follows: 

5°  = 25, 7°  = 4,  R°  = 6,  r = .4,  J3  = 0.2,  a = 0.09,  p = 0.3,  p = 0.3,  S = 0.2,  t = 1,  K = 50. 

Then  we  get  R0  = 1.136  >1,5  =11,7  = 1.21,  by  theorem  3.3, 5(f)  and  7(f)  approach  to  their  steady-state 
values  with  delay,  the  disease  will  be  exist. 

5.  CONCLUSIONS 

We  have  analytically  studied  a delayed  SIR  model  with  the  exponential  demographic  structure  and  the  nonlinear 
incidence.  We  found  the  sufficient  condition  of  the  stability  for  the  endemic  and  disease  free  equilibrium  of  the  model. 

When  R0  < 1,  the  disease  free  steady  state  is  stable,  no  other  equilibria  exist.  When  R0  > 1,  a unique  endemic  equilibrium 
exists  and  stable  under  some  condition  with  and  without  delay. 
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